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1. Introduction 

Lattice study of the p meson decay is an important step for understanding of the dynamical 
aspect of hadron reactions induced by strong interactions. There are already three studies |JIJ ^J, ^J. 
The earlier two studies employed the quenched approximation ignores the decay into two ghost 
pions. The most recent work, while using the Nf = 2 dynamical configurations, concentrated on 
the p — ► KK transition amplitude rather than the full matrix including the p — > p and %% — > Tin 
amplitudes. All three studies were carried out at an unphysical kinematics m K /m p > 1 /2. 

In this work we attempt to carry out a more rigorous approach. We estimate the decay width 
from the scattering phase shift for the 1 = 1 two-pion system. The finite size formula presented by 
Rummukainen and Gottlieb [Q] is employed for an estimation of the phase shift. Calculations are 
carried out with Nf = 2 full QCD configuration previously generated for a study of light hadron 
spectrum with a renormalization group improved gauge action and a clover fermion action at j8 = 
1.8, K = 0.14705 on a 12 x 24 lattice [Q]. The parameters determined from the spectrum analysis 
are l/a = 0.92 GeV, m K /m p = 0.41, and L = 2.53 fm. All calculations of this work are carried 
out on VPP5000/80 at the Academic Computing and Communications Center of University of 
Tsukuba. 



2. Method 

In order to realize a kinematics such that the energy of the two-pion state is close to the 
resonance energy m p , we consider the non-zero total momentum frame (the moving frame) [Q] 
with the total momentum p = In/L-e^. The initial p meson is assigned a polarization vector 
parallel to p. One of the final two pions carry the momentum p, while the other pion is at rest. The 
energies ignoring hadron interactions are then given by = \J m| + p 2 + m % for the two-pion state 



and Z?2 = y m p + P 2 f° r tne P meson - We neglect higher energy states whose energies are much 
higher than E® and E®. On our full QCD configurations, the invariant mass for the two-pion state 
takes yfs = 0.97 x m p , while 1.47 x m p is expected for the zero total momentum. The p meson at 
zero momentum cannot decay energetically, so that it can be used to extract m p . 

The hadron interactions shift the energy from Zs„ to E n (n = 1,2). These energies E n are re- 
lated to the two-pion scattering phase shift 8(y/s) through the Rummukainen-Gottlieb formula [Q], 
which is an extension of the Luscher formula ^ to the moving frame. The formula for the total 
momentum p = pe^ and the representation of the rotation group on the lattice reads 

1 1 l + {3r 2_ r 2 )/q 2 

tanS(V^) 27i 2 qy £f r 2 -q 2 ' K ' ' 

where y/s = \/E 2 — p 2 is the invariant mass, k is the scattering momentum (y/s = 2\/m 2 +k 2 ), y 



is the Lorentz boost factor (7 = E/y/s), and q = kL/ (2n). The summation for r in (2.1 ) runs over 
the set 

T = {r| r 1 =n 1 ,r 2 = n 2 ,r 3 = (n 3 + ~)/y , n G Z 3 } . (2.2) 

2 27t 

The right hand side of ( |2.1[ ) can be evaluated by the method described in Ref. [f7|]. 



2 



p meson decay from the lattice 



N. Ishizuka 



In order to calculate E\ and E2 we construct a 2 x 2 mauix time correlation function, 



= ( (o\ (*»)(<.) |o) (o| (^)+(0 P3 fe) |o) 

U I (0| p 3 t (f) (nn)(t s ) |0) (0| />}(*) p 3 (4) |0> 



(2.3) 



Here, pi,(t) is an interpolating operator for the neutral p meson with the momentum p = 2n/L ■ e 3 
and the polarization vector parallel to p; (nn){t) is an interpolating operator for the two pions given 
by 

(K7l)(t) = ^=^-(p,t)7l + (0,t)-7l + (p,t)7l-(0,t)') , (2.4) 

which belongs to the and iso-spin representation with / = 1, I z = 0. 

We can extract the two energy eigenvalues by a single exponential fitting of the two eigenvalues 
Ai(f,f#) and A2 (? , ) of the normalized matrix M(t,t R ) = G(t)G~ l (f#) with some reference time 
tR [§] assuming that the lower two states dominate the correlation function. 

In order to construct the meson state with non-zero momentum we introduce a U(l) noise 
§/(x) in three-dimensional space whose property is 



1 N R 

-£y(x)^(y) = 5 3 (x-y) for 



We calculate the quark propagator 

Q(x,t\q,t t ,Zj) = ^(D- 1 ) (x,/;y, t s ) 



Nr 



(2.5) 



(2.6) 



regarding the term in the square bracket as the source. The two point function of the meson with 
the spin content Y and the momentum p can be constructed from Q as 



1 N * 1 



(2.7) 



where the bracket refers to the trace for color and spin indeces. 

The quark contraction for the %% — > %% and the %% — > p components of G(t) are given by 



Gtyty— >7T7T (^) 



-P 



P 

-p 





P 



+ 




+ 





(2.8) 

where the four verteces for the %% — > 7T7T and three verteces for the 7T7T — > p components refer to 
the pion or the p meson with definite momentum. The time direction is upward in the diagrams, 
and the p — > 7T7T component is given by changing the time direction. 
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The first term of the Tin — ► %% component in (ETHJ) can be calculated by introducing another 



U(l) noise rjj(x) having the same property as £y(x) in ( |2.5| ); 



R 7=1 *.y 



^L£e-*- x -((2 + (x,f|<U,^) 0(x ) f|p,f„^)>(0 t (y»*|O,f, ) ^) 0M<M*>ty)) . (2.9) 



The second term of ( |2.8| ) is obtained by exchanging the momentum of the sink in ( |2.9| ). 



In order to construct the other terms of ( |2.8| ) we calculate a quark propagator of another type 
by the source method, 



W(x,t\k,t 1 \q,t s ,^j)=Y J (D l )(x,t;z,h) ■ ie**)* gfofilq,^) 



(2.10) 



where the term in the square bracket is regarded as the source in solving the propagator. Using W 



we can construct the third to sixth terms in the 7171 — > 7T7C component of ( |2.8| ) by 
3rd 



1 N R 



4th 
5th 
6th 



w\x,t\0,t s \ -p,t s ,£j) W(x,t\0,t\0,t s ,£j 
W(x,t\0,t s \p,t s ^j)W\x,t\0,t\0,t s ,^j 



1 N s 

1 Nr i 

— £JVW W (^ t \ pts \o t ^) W ^ x ^ t \0,t s ^j 



1 



Nr 



y v e -* x 



W^rl -p,f,|0,f„&) W(x,f|0,f|0,f„£y) 



(2.11) 



The two terms of 7C7C — > p of ( |2.8[ ) can be similarly constructed by 
1st 



Nr 



2nd 



1 N R 



ipx 



W f (x,t\ -p,t s \0,t s ,Zj) (75B) Q(x,t\0,t t ,£j) 



E£e-*- x -( fit^fio,,,,^) ( 7573 )W(x,t\p,t s \0,t s ^j) 



(2.12) 



In this work we set the source at t s = 4 and impose the Dirichlet boundary condition in the 



time direction. We calculate the <2-type propagators for four sets of q and the U(l) noise in (2.6) : 
(q, noise) = {(0,^),(0,T]),(p,^),(— p,£)}. The VK-type propagators are calculated for 22 sets 
of k, t x and q in : (k,/i|q) = {(p,f,|0), (-p,f,|0), (0,f,|p), (0,f,| -p), (0,^ = 4 - 21 10)}, 

with the same U(l) noise All diagrams for the time correlation function can be calculated with 
combinations of these propagators. We choose Nr = 10 for the number of U(l) noise. We carry out 
additional measurements to reduce statistical errors using the source operator is located at t s + T/2 
and the Dirichlet boundary condition is imposed at T/2. We average over the two measurements 
for the analysis. Thus we calculate 520 quark propagators for each configuration. The total number 
of configurations analyzed are 800 separated by 5 trajectories 

3. Results 

In Fig. [I] we plot the real part of the diagonal components (7T7r — ► 7T7T and p — > p) and the 
imaginary part of the off-diagonal components {%% —> p, p — » 7T7T) of G(t). Our construction of 
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Figure 1: G(t) 
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Figure 2: Normalized eigenvalues X\(t,t R ) and X2(t,t R ). 
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G(f) is such that the sink and source operators are identical for a sufficiently large number of the 
t/(l) noise. In this case we can prove that G{t) is an Hermitian matrix and the off-diagonal parts 
are pure imaginary from P and CP symmetry. We find that this is valid within statistics, but the 
statistical errors of the p —> %% component is larger than those of %% — > p in Fig. p]. In the following 
analysis we substitute p — > %% by %% — > p to reduce the statistical error. 

The two eigenvalues Ai (?,?#) and A2(/,/r) for the matrixM(f,f«) = G{t)G~ x {Ir) are shown in 
Fig. ||. We set the reference time ?r = 9 and normalize the eigenvalues by the correlation function 
for the free two-pion system, (0|tt(— p,f)w(p,f,)|0) (0|tt(0, ?)tt(0, ^ ) |0). Thus the slope of the figure 
corresponds to the energy difference AE„ = E n — E®. We observe that the energy difference for X\ 
is negative and that for A2 is positive. This means that the two-pion scatting phase shift is positive 
for the lowest state and negative for the next higher state. 

We extract the energy difference AE n for both states by a single exponential fitting of the nor- 
malized eigenvalues X\ and A2 for the time range t = 10 — 16. Then we reconstruct the energy E„ 
in the moving frame by adding the energy of the two free pions, i.e., E„ = AE n +E®, and convert it 
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Figure 3: sin 2 5 (•>/$), position of m p and resonance mass Mr. 



to the invariant mass y/s. Substituting y/s into the Rummukainen-Gottlieb formula (2.1) we obtain 
the scattering phase shift : 



a\f~s 


tan S(\/s) 


0.7880 ±0.0082 


0.0773 ±0.0033 


0.962 ±0.024 


-0.43 ±0.12 



(3.1) 



The p meson mass obtained at zero momentum is am p = 0.858 ±0.012. Hence the sign of the 
scattering phase shifts at y/s < m p is positive (attractive interaction) and that at y/s > m p is negative 
(repulsive interaction) as expected. The corresponding results for sin 2 8 (y/s), which is proportional 
to the scattering cross section of the two-pion system, are plotted in Fig. || together with the position 
of m p . 

In order to estimate the p meson decay width at the physical quark mass we parameterize the 
scattering phase shift by the effective p — > nn coupling constant g pnn , 



tanS(^) = ^ 



k 3 



6k yfs{M 2 R -s)' 
with g pK7[ defined by the effective Lagrangian, 

L e s. =g p wE abc (k l -k 2 ) ll p^{p)n b {ki)n c {k2) , 



(3.2) 



(3.3) 



where k is the scattering momentum and Mr is the resonance mass. The coupling g p7[K generally 
depends on the quark mass and the energy, but our present data at a single quark mass do not 
provide this information. Here we assume that these dependences are small and try to estimate 



gpmt and Mr from our results in (3.1 ). We also estimate the p meson decay width at the physical 
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quark mass from 

r " = ! 6? ■^=£™ x4 - m MeV < ( 3 - 4 ) 

where m p is the p meson mass at the physical quark mass and k p is the scattering momentum at 
Vs = rhp. 

Our final results are as follows. 

aM R = 0.877 ±0.025 
g p7lK = 6.01 ±0.63 
T p = 149 ±31 MeV . (3.5) 

The resonance mass Mr obtained from the scattering phase shift is consistent with am p = 0.858 ± 
0.012 obtained from the p meson with zero momentum. The p meson decay width T p at the 
physical quark mass is consistent with experiment (150 MeV). In Fig. || we indicate the position of 



Mr and draw the line given by (3.2) with g pK7[ and Mr in (3.5) 



4. Summary 

We have shown that a direct calculation of the p meson decay width from the scattering phase 
shift for the 1 = 1 two-pion system is possible with present computing resources. However, several 
issues remain which should be investigated in future work. The most important issue is a proper 
evaluation of the quark mass and energy dependence of the effective p — > nn coupling constant 
gpxK- This constant is used to obtain the physical decay width at m K /m p =0.18 from our results at 
m n jm p = 0.41 by a long chiral extrapolation. In principle we can estimate the decay width from 
the scattering phase shift without such a parameterization, if we have data for several energy values 
at or near the physical quark mass. This will be our goal toward the lattice determination of the p 
meson decay. 

This work is supported in part by Grants-in-Aid of the Ministry of Education (Nos. 
17340066, 16540228, 18104005, 17540259, 13135216, 18540250, 15540251, 13135204, 
16740147, 18740139 ). The numerical calculations have been carried out on VPP5000/80 at 
Academic Computing and Communications Center of University of Tsukuba. 
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